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Abstract 

Using finite gap methods, we find the leading order finite size corrections for 
an arbitrary number of giant magnons on physical strings, where the sum of the 
momenta is a multiple of 27r. Our results are valid for the Hofman-Maldacena 
fundamental giant magnons as well as their dyonic generalizations. The energy 
corrections turn out to be surprisingly simple, especially if all the magnons are 
fundamental, and at leading order are independent of the magnon flavors. We also 
show how to use the Bethe ansatz to find finite size corrections for dyonic giant 
magnons with large iZ-charges. 
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1 Introduction 

Integrability in A/" = 4 planar gauge theories [1,2] has provided many deep insights into 
the AdS/CFT correspondence [3]. Recent spectacular progress [4,5] in finding expressions 
for infinite dimensional operators valid for all values of the coupling rely on the existence 
of an integrable S-matrix [6-11] for a spin chain with infinite extent and long range 
interactions. 

Still, there remain many unsolved problems. Chief among these is the question of 
what happens for operators of finite size. In this case, wrapping effects arise due to 
the long range nature of the interactions [12], and appear to spoil the Bethe ansatz 
equations in [7]. Nonetheless, there has been some progress in this direction. In [13], 
Ambjorn, Janik and Kristjansen gave a systematic analyis of wrapping effects using the 
Thermodynamic Bethe ansatz. Their analysis showed that wrapping effects will first 
appear at the L loop level for an operator which corresponds to a spin chain with L 
sites, precisely as expected. Further analysis was carried out by Kotikov et. al. [14] who 
argued that wrapping effects must be included, lest there be a breakdown in the BFKL 
constraints at 4 loop order. There have also been some recent explicit computations 
of the four loop contribution to the Konishi operator anomalous dimension [15, 16], 
where wrapping effects lead to a C(5) term, a factor not present in the asymptotic Bethe 
equations 

On the string theory side, finite size effects were computed for spinning strings by 
Schafer-Nameki, Zamaklar and Zarembo [17]. Here they considered a large spin S and R- 
charge J and found that the asymptotic Bethe equations in [7], including the appropriate 

^Whilc both [15] and [16] find a C(5) contribution, there is a discrepancy between the two results. 
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dressing factors in [10] and [11] give the correct answer for the energy levels up to terms 
of order exp(— 27r J/v^A), where A is the 't Hooft coupling. 

Another direction for investigating finite size effects is to compute the corrections for 
giant magnon [18] energies. On a chain where the i?-charge J is taken to infinity, the 
energy of a single fundamental magnon on the resulting infinite chain is [8] 



£= ^l + lQgHm^l, (1.1) 
where g is related to the 't Hooft coupling by 



and p is the world-sheet momentum of the magnon. In the classical world-sheet limit 
where g ^ 1, the corresponding giant magnon solution was found by Hoffman and 
Maldacena, where the classical result is 

^^4^ sin 2. (1.3) 
One can also consider Q bound magnons [19], where the resulting energy is 



8 = jQ^ + lQg^sin^l. (1.4) 



The quantity Q is the value of a second i?-charge, and if Q ^ 1, the result in (1.4) can be 
found classically, the so-called dyonic giant magnon [20]. One can think of the classical 
result in (1.3) as the limiting case of (|L4| when (5 = 00 



One can then consider corrections to the energies of ( |1.3[ ) and (1.4) for finite but very 



large J, the criteria being that J ^ g and J ^ Q. Such a calculation was carried 
out by Arutyunov, Frolov and Zamaklar, using a light cone gauge-fixing on the string 
world-sheet, where they found that the leading correction in the (5 = case is [22] 

Z\£ ~ -16^ sin^fe-^exp ( -2 — ^—-\ (1.5) 



(An implicit result was also given in [22] for the (5 7^ dyonic case.) The result in (1.5) 



is for a single magnon. In an entirely different approach, Janik and Lukowski derived 



the correction in (1.5) using Luscher's method for finite size effects in a relativistic 
quantum field theory. The approach of these authors required the full BHL/BES S- 
matrix, providing a further check of the dressing phase in [5] . 

However, strictly speaking there cannot be only one magnon, since in this case the 
total momenta of the magnons is not a multiple of 27r, which violates the Virasoro 
constraints. From the point of view of the dual gauge theory, having a total momenta 
that is not a multiple of 2tt violates the trace condition. Hence, for a single magnon one 
would expect gauge dependence in the final answer. Indeed in [22] it was shown that the 



= Q is the classical charge of a fundamental giant magnon. The actual charge is Q = 1, which is 
a consequence of fermion zeromodes [18,21]. 
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higher order correction terms were dependent on a gauge parameter. One way around 
the predicament is to put the theory on a Z„ orbifold. In this case the total momentum 
only needs to be a multiple of li^jn. This was done in [23] where they confirmed the 



result in (1.5). Or one could just have n copies of the magnons, so that one finds the 
limit of the GKP spinning string on [24] or its multicusp generalization [25]. 

But another way to proceed is to simply add more magnons, such that the total 
momentum is mod 27r. With the Virasoro constraints now satisfied the results should be 
completely gauge independent. However, in this case the presence of the other magnons 
would be expected to modify the correction to each magnon's energy. 

In this paper we will explicitly compute the leading order corrections to the magnon 
energy in the presence of other magnons. We also include the possibility of bound 
magnons in the analysis. Our main method employs the finite gap method in [26] for the 
string cr-model constrained to an i? x S*^ subspace. S'^ has an SU{2) x SU{2) isometry 
and so our classical string configurations come with two independent i?-charges. 

Solutions of the finite gap equations will have 2M square root branch cuts, where M 
is the number of giant magnons (some or all of these giant magnons could be dyonic). 
Hence the general solutions will involve genus 2M — 1 hyperelliptic functions. However, 
in the limit that J — oo the hyperelliptic surface is highly singular with the branch cuts 
shrinking off to zero size. In this limit it is possible to make a vast simplification and 
approximate the hyperelliptic functions with ordinary trigonometric functions. One of 
our main results is that for M fundamental giant magnons, the energy correction for 
each magnon is 

M ■ 2 Pj+Pk 

ASj = -16 g sin^ f J] ^^^T^ exp(-2( J + ^ S,)/Sj) < p, < 2n (1.6) 

kf^j 4 k 

Note that the expression J + ^ is the total energy E to leading order and can be 
thought of as the natural length of the chain since each magnon has a size equal to its 
energy. Hence the natural suppression factor is e~^^/^J , not e"^'^/^^ . 

We can also extend our results to nonzero values of Q, although the final expressions 



will not have quite the elegant form as in (1.6). Moreover, the exponential suppression 
term is modified to 

where 

g^sin^f + 16(72sin4| 

It is now also possible to compare these results directly with the Bethe ansatz. In the 
limit where all Qj ^ g, the string result will approach that found for the one-loop 
gauge theory results. In this case the problem reduces to finding the spectrum of the 
Heisenberg spin chain, where one can explicitly apply Bethe's ansatz. The computation 
is significantly simpler but will still agree with the finite gap result. 

In section 2 we show how to derive finite gap equations from the Heisenberg spin- 
chain and then give the modifications for the Rx cr-model, referring the reader to the 
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references for the details of this latter part. We then present the general configuration for 
a collection of giant magnons at large but finite J. In section 3 we solve the gap equations 
to leading order for a generic giant magnon configuration. We then go through some of 
the specific cases. In section 4 we show how to find similar results for the Heisenberg 
spin chain and how these match onto the solutions in the previous sections. In section 
5 we discuss how to reproduce the results for a single set of bound magnons using TBA 
analysis, obtaining the same exponential term. In section 6 we present our conclusions 
and suggestions for further work. 

The results appearing here were first presented in [27]. As this paper was being 
typed, we received [28] which considers the finite size corrections for a single dyonic 
giant magnon. 



2 From Heisenberg to finite gap 

In this section we give a short review for applying the Bethe ansatz to the one-loop 
anomalous dimension for single trace operators with large i?-charge. We show how to 
derive a set of integral equations and then generalize these equations to the finite gap 
solutions of the string sigma model on R x . 

The starting point is the chiral primary operator Tr Z"^ , whose dimension is protected 
by supersymmetry. The field Z is a complex adjoint scalar and the it!-charge J is a 
conserved global charge. We can build other operators by inserting impurities between 
the Z-fields to give 

Tt[Z . . . ZxiZ . . . Zx2Z . . . ZxsZ ...Z] (2.1) 

where the Xi ^-^e one of eight boson or eight fermion fields. For what follows we will 
restrict ourselves to the SU (2) sector, so the operators we consider will be of the form 

Tr[Z . . . ZWZ . . . ZWZ . . . ZWWZ . . . ] 
- Tr[T ... Tit ■■■ Tit ■■■ Tiit ...], (2.2) 

where we have shown the map between operators of this type and a spin chain with 
up and down spins. In [1] it was argued that the one-loop anomalous dimensions for 
operators of this type are equivalent to the energies of the Heisenberg spin chain with 
Hamiltonian 

H^L + 2g^J2(^l-'^^i- -^^+1) (2-3) 

where L is the total number of sites on the lattice. If we assume that there are J Z-fields 
and Q W^-fields then L = J + Q. Each of the i?-charges J and Q are conserved, and will 
be conserved to all orders in the perturbative expansion. Note that we have included a 
constant term L in the Hamiltonian so that the energies will match the full dimension 
up to the one-loop level. 

The Heisenberg spin-chain is an integrable system and can be solved by finding a set 
of momenta for the impurities on the chain. The impurities are the down spins, in other 
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words the Ty-fields. Each momentum satisfies a quantization condition 



Q 

e^P^^'+^) = \[S-\p,,pu) (2.4) 

where S{pj,pk) is the S-matrix for impurities j and k and we write the number of lattice 
sites as a sum of the number of Z and W fields for later convenience. It is convenient to 
express the momentum in terms of a rapidity variable Uj, where 

■ (2.5) 



in which case the quantization condition in (2.5) becomes the well known Bethe ansatz 
equations 

I • /o\ Q+J Q I • 

Uj-i/2) l}^Uj-Uk-i' 

The Uj are sometimes referred to as "Bethe roots". The energy of a particular state is 
additive with the individual energies of the magnons and has the form 

Q 

i=i 

where 

The momenta pj must also satisfy one further condition that one does not normally 
encounter for the spin chain. Since the gauge invariant operator has a trace, it is invariant 
under a single shift of all the scalar fields inside the trace. This translates into the trace 
condition for the magnon momenta 

Q 

"^Pj = mod 27T . (2.9) 

i=i 

Now let us go into the thermodynamic limit and assume J ^ Q ^ 1. If we take the 



log on both sides of the Bethe equations in (2.6) we have the approximate equation 



.Q + J 

I 



^3 t^j ^3 - «fc 



where the integer rij arises from choosing a branch of the log. If we define the density 

p(x) = ^(5(x-Mfc) (2.11) 



k 
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then the equation in (2.10) becomes an integral equation [1,29] 



dx'p{x') Q + J 



2+ — '-^-^ = 27rn,- 



(2.12) 



where the contour Cj refers to the contour where all roots are on the same log branch 
rij. Clearly the density must satisfy 



dx'p{x') = Q 



{J + Q)-J + Q 



(2.13) 



where the second expression is given for later comparison with the cr-model. 

In the thermodynamic limit we can approximate the momenta of the individual magnons 

as 



Pj - V 



Hence, it follows that the total momentum P is given by 

dx' p{x') 



x' 



P = mod 2n 



(2.14) 



(2.15) 



where the last equality follows from the trace condition. Furthermore the energy of each 
magnon may be approximated as 



~ 1 + 2^V«? , 

from which one finds the total energy 

dx'p{x') 



E = J + Q + 2g^ 
We will rewrite this equation as 



X 



dx'p{x') _ 1 E- J -Q 



X 



9^ 



(2.16) 



(2.17) 



(2.18) 



Hence, our task is to find solutions for p(x) in the integral equations in (2.12), while 
satisfying the conditions in (2.13) and (2.15). The integral equation was found by taking 
the classical limit J,Q ^ oo. As such, the solutions should correspond to solutions of 
some classical system. Indeed the classical system is the Landau-Lifschitz model which 
has the equation of motion 

:S = -2g^Sx—S, (2.19) 



dt 



where a is the one dimensional position. The spin S should be thought of as a continuous 
variable with fixed length |. 

Let us now turn to the general solutions of (2.12). We start by defining the resolvent 

Gix), 



G{x) 



dx'p{x') 
X — x' 



(2.20) 
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Figure 1: A solution with 4 branch cuts. The different cuts correspond to different log branches 
and are symmetric under complex conjugation. 



where the integral is over N contours where p{x') has support and x is assumed to be 
off of the contours. We will assume that p{x') is zero at the end points of the contours 
and so the resolvent will have square root branch points at these endpoints, with the 
branch cut giving a discontinuity across the contour. A typical example of the cuts is 
shown in figure 1. Note that the cuts are symmetric under complex conjugation which 
is a consequence of the Bethe equations and guarantees real values for the total energy. 

We can then extend the resolvent G{x) to a two-sheeted surface connected by the 
branch cuts. Hence this surface has genus — 1 and is hyperelliptic. The top sheet is 



called the "physical" sheet and from (2.20) it is clear that G{x) is nonsingular and single 



valued on this sheet. Using the integral equation (2.12) we see that along a contour 
G{xj) satisfies the condition 

G{xj + ie) + G{xj - ie) = 9J:I. _ 2nnj xj G Cj (2.21) 

Xj 

Hence, we should expect G{x) to have a pole on the bottom sheet at x = 0. We can also 
define two types of cycles, A- and B-cycles, where A-cycles circle around single cuts on 
the top sheet and B-cycles cross one cut onto the bottom sheet and cross another cut 
back onto the top sheet. Such cycles are shown in figure 2. If we integrate dG around 
an A-cycle, we find 

j dG = Q (2.22) 

because of the single-valuedness of G. On the other hand, integrating dG around a 
B-cycle that traverses cuts j and k gives 

dG = 2n{nj - Uk) (2.23) 

B 
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A cycle 



Figure 2: An A-cycle circles around a single branch cut, while a B-cycle traverses two cuts in 
going from the top sheet to the bottom sheet and back again. 



which follows from (2.21). 



Now it is possible to recut the surface such that the hyper surf ace is unchanged. This 
leaves the branch points where they are but changes how the branch points are paired up 
onto branch cuts. These transformations correspond to Sp{N — 1,Z) transformations of 
the hypersurface. These transformations should leave the integrals around the A and B- 
cycles unchanged, so this forces us to introduce condensates. A condensate extends along 
the imaginary direction with a constant density which is quantized, namely p{x) = —in 
where n is a positive integer. The factor of —i is so that / dxp{x) is positive real. In figure 
3 we recut the surface in figure 1 so that two branch cuts are redrawn. A condensate now 
extends between the upper and lower branch cuts so that the B-cycle remains nontrivial. 



If the B-cycle satisfies (2.23) then the condensate density is 

p(x) = -i{nj - rik) . (2.24) 

In order that the original A-cycles around the j and k cuts still give zero, we further 
require that the resolvent satisfy 

G(x + ie) + G(x - le) = 9AZ - 2nn xeC (2.25) 

X 

where C is either of the new contours and n is some integer. Note that n has to be the 
same on both contours so that the A-cycle integrals are zero. There is some freedom to 
move the condensate along the new branch cuts, but for simplicity we will place its end 
points on two of the branch points. In order that there be no new poles in the resolvent, 
it is necessary that the contour density at these branch points equals the condensate 
density. 

Suppose for now that there are two contours and recut the surface so that there is a 
condensate. Let us further assume that the density of the condensate is p{x) = —i. In 
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Figure 3: Recutting the hyperelliptic surface. The condensate is now added to preserve the 
cycles. 



general this configuration will not satisfy the trace condition, but this is certainly not 
a problem for the Heisenberg spin chain. The resulting contours are shown in figure |4] 
where the branch points are labeled by A, B and their complex conjugates. 

Let us now take the limit J ^ oo. In this limit the branch points A and B approach 
each other and the surface degenerates [30]. In this limit we are left only with the 
condensate [31,30]. It is then straightforward to show that 

Q = j dxp{x) = -i{A- A*) 

p = /^ = -.iog(AM*) 

J X 

E-J-Q = 2g' [^^^ = +2ig^A-'-A*-'), (2.26) 



from which it follows 



A = |e^P/2cscf 



,2 



£ = E-J = Q + 8yr sin^ 2 ^ Jq2 + 16^2 ^^^2 p_ (2.27) 

where we have included the approximate expression ioi E — J for comparison with the 
(T-model. The expression £ corresponds to the energy of Q bound magnons with total 
momentum p. Note that if we had chosen the condensate density to be p{x) = —i n, 
then the resulting solution is n sets of bound magnons, each with charge Q/n and each 
with momentum p 



The classical solution that (2.26) corresponds to is the soliton wave discussed by Lak- 



shmanan, Takhtajan and Fogedby in the classical one-dimensional ferromagnet [32]. The 
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X 



Figure 4: Two cuts joined by a condensate with density p{x) 



OB 



p(^)=-i 



soliton is a localized spin wave with a width and velocity dependent on the momentum. 
The polar direction of the spin as a function of time and position a is given by 



6 = 2 sin"^ 
where the width F and the velocity v are 

Q 



sin I sech — 



r 



2sin2 2 



Q 



smp 



(2.28) 



(2.29) 



Let us now generalize the previous results for the Heisenberg spin chain to that of the a- 
model on i? X 5*^ with Virasoro constraints. The details of how the integral equations are 
derived can be found in [26], so here we only present the results. First, the generalization 



of the integral equation in (2.12) is 



dx'p{x') 



Ex 



x^ 
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- - 27rnj 



X e Cj 



(2.30) 



where we see that the difference between (2.12) and (2.30) is that the pole at x = has 



split into two poles at x = ±g, and the residue has changed from J + Q to E. The only 



other change we need to make is that (2.13) is modified to 

dx' p{x') = ~^ Q ^ 



(2.31) 



where like in (2.30), the difference between (2.30) and (2.12) is that J + Q is replaced 



with E. The expressions (2.15) and (2.18) are unchanged 
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If we now take the two cut solution, recut with a condensate, and then take J ^ oo, 
we are left only with the condensate [31,30]. In this case we find 

^{E-J + Q) 
P 

E-J-Q = +2ig\A-^ - A*~^) , (2.32) 



-i{A-A*) 

-i\og{A/A*) 

+ 2ig^ {A-^-A*~^] 



from which we are lead to 



S = E-J 



g2 + 16^2 sin2 2 



(2.33) 



The classical solution that this condensate corresponds to is the Hofman-Maldacena 
giant magnon [18] if Q = or its generalization if Q 7^ [20] . Note that in the limit that 
g ^ (2.33) reduces to (2.27) for the bound magnons in the Heisenberg model. 



3 Finite size corrections for giant magnons 



In this section we find the leading order finite size corrections for the energy of a giant 
magnon in the presence of an arbitrary number of other giant magnons. With more than 
one magnon it is possible to satisfy the Virasoro constraint that the total momentum is 
mod 2n, so the results will be gauge independent. 

The root configuration corresponding to several finite size giant magnons is depicted 
in figure |5| A number of cuts are pairwise connected by condensates. To calculate the 
dispersion relation, we want to solve the equations (2.15), (2.18), (2.31) and (2.30). Even 



for a single magnon, the solution of (2.30), which leads to the AFZ result, is given by an 



elliptic function. For several magnons hyperelliptic functions are needed. The problem 
is, however, significantly simplified if we assume each cut to be very short compared to 
the condensates. 



For X on a given cut AB, we can write the integral equation (2.30) as 



^ dx'p(x') x-A r^* dx'pix') 

T + 2 « In — + / 

X — x' X — A* Jq, 



X — x' 



dx'p{x') 

others 



E 



+ — 1 - 2?™. (3.1) 



2 \x — g X + g 

Except for the first two terms on the left hand side, this expression depends only weakly 
on X. Thus, to a reasonable approximation. 



B 



dx' p{x') 



-2i ln{x-A) + C, 



f X — 

where C is a constant. By evaluation a.t x = A, C can be expressed as 

C ^ — I -r^ + -^—] -2Tm + 2i \n(A- A*) -2i Vln^""^' 



A- A* 



(3.2) 



(3.3) 
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Figure 5: Cuts for several finite size giant magnons 



The problem has now been reduced to an integral equation over a single contour, which 
can be solved for p using a finite inverse Hilbert transform [33], 

_ ^iB-x)ix-A) dx' -2t\n{x'-A) + C 
2 h x-x' ^{^B-x'){x' -A) 

The term that involves the constant will vanish, while the rest of the integral can be 
calculated to give 

/ N 2i IB- A , , 

pix) = arccos \ -. (3.5) 

TT \ X — A 

Note that p{B) = while p{A) = —i. Hence p{x) is continuous where the cut joins the 
condensate, so the resolvent (2.20) will not have an extra pole there. Such poles should 
only occur in the limit A where the surface degenerates [30] |^ 

In the J — > oo limit, the end point of the condensate is at 

A Q + S iE P , . 

Ao = xo + lao = — - — e 2 esc - (3.6) 

For finite J this end point is shifted to A = Ao + 6 A = xo + ia and the end point of the 
cut is ai B = xq + ib. Since the cut is assumed to be small we write 

B-A = i{b-a) = i6e'''', (3.7) 



with 5 -C a. Plugging p{x) back into (3.2), we can now perform the integral to calculate 
C, which will be of order 0{S). Comparing this to ( |3.3[ ) we express the shift B — A as 

Se^'^ = 8e-^(^+^)+*-'^n 4^- (3-8) 

j ■' 



*We thank N. Dorey and B. Vicedo for discussions on this point. 
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For J — > oo, the charges p and Q are given by 

p = -z In g = -^(^ - ^o) - - A*-'). (3.9) 

For finite J there are additional contributions from the change of the integral over the 
condensate, due to the shift of the end point, and from integrals along the cuts. Requiring 
that the shifts of p and Q vanish to the leading order, leads to the equations 

, 5A 5A*\ 
-M^-^ -0, (3.10) 













e~ 




^ A^ 



A*oJ \Ao Al 



2a /42_«2 4*2 _ «2 



'^cos^ - ^ ( ^ + ^ ) - . ( ^^5A - ^^5A* ) = 0, (3.11) 



r, \ a2 ' 4*2 \ A2 a* 

^ ^0 / \ ^0 ^0 

which is solved by 

6A = ~^-6e''^. (3.12) 

We see that the end point of the condensate splits into a cut. The shift in the dispersion 
relation is given by 



The first nonvanishing correction occurs at order 5^. This correction is calculated in the 
same way — the contributions to the charges from the shifted contours are calculated 
and required to vanish. The resulting expressions for the shift of A is given by the 
somewhat complicated expression 



2 S{S + QY 



Sg' + {£ + Qf . 

^- — sm(p — 20) 



2sin(3p-2(/)) -4sin2(/) + sin(p + 20) ^ • 4 P 

g — ^ ^ + Aig cos{p — 20) sm - 



+ 0{5^). (3.14) 



The quadratic order correction to the dispersion relation, however, simplifies to 



2 

A£ = --^ sin^ ^ cos(p - 20) 5\ (3. 15) 



We will now consider this expression in several cases. 
3.1 Single fundamental magnon (Q = 0) 

First we consider a single fundamental magnon with momentum p. In general such a 
configuration will not satisfy the momentum condition. The expansion parameters 5 and 
are given by 

5e''l' = 8aoe-^'^"exp(-E/£) = Saoe"'™ exp(-( J + £)/£). (3.16) 
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Insertion into (3.15) gives 



AS 



■ ^ sin^ ^ cos(j9 + 27m) exp(— 2J/£^) 
c 2 



sin 



3 V 



cos(p + 27m) exp I —2 



J 



Ag sin | 



(3.17) 
(3.18) 



Comparing this to the AFZ result we see that, for integer n, the results disagree by 
a factor of cos p. Of course, if we allow n to be noninteger then the cosp factor can be 
removed. In fact this noninteger n is a remnant of the missing momentum necessary to 
satisfy the level matching condition. We will see shortly that with additional magnons, 
their momenta will enter into the argument of the cosine. 



3.2 Single dyonic magnon {Q ^ 0) 

For non-zero Q the parameter 6 and the phase are changed to 

Ocot 2 

S = 8aQ exp{-ES/S^), (j) = E - tth, (3.19) 

^; 

with ^ 

S^^-^ + lGg'sin'^. (3.20) 
sm I 2 

This results in a shift of the dispersion relation of the form 

2 

AS = -64sin^| cos(p-20)-|^exp(-2J£:/£,^), (3.21) 

where 

"-^2- g2 + i6^?2 sin^ 2 • ^^■^^) 

Note that for E ^ Ss this correction is highly oscillatory. 

Further note that we can get rid of the e~^" term if we replace J by E. The interpreta- 
tion of this is clear. The size of a magnon is basically its energy, so we would expect the 
total size to be J + £^ = -E. Hence the exponential suppression factor should be governed 
by the size E and not J. 

3.3 Several fundamental magnons (all with Q = 0) 



Equation (3.8) now reads 

S, e^^^ = 8ao, n exp{-E/£,) , = Ag sin |, (3.23) 

k^j ^ 
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where our notation is extended in an obvious way to allow for several parameters S, 0, 
ao and p. Using the lowest order solution for Ak and assuming < pa; < 27r we have 



exp 




sm 



sm 



(3.24) 



resulting in 



AS, 



AQgsin^^Yl 



M - 2 Pj+Pk 
3 Pj TT 4 



2 Pj -Pfc 



sm 



exp(-2(J + ^^fc)/£:, 



(3.25) 



Let us note a few points about the result in (3.25). 



For finite j the momenta will be quantized and this quantization will depend on the 
magnon S-matrix. Hence (3.25) applies for the pj that satisfy this quantization. 



Magnons come in flavors, but the result in (3.25) is flavor independent at this level. 
In fact the flavors arise because of fermion zeromodes [21], which is obviously a 
quantum effect so one should not expect to see their effect in a classical result. In a 
fully quantized theory we would therefore expect the flavor effects to be subleading 
in 1/g. 



The result in (3.25) is singular if pj pt- In this limit we are left with two 
magnons with identical momenta. We will see the resolution of this in the next 
subsection. 



3.4 M magnons with the same momentum and Q 

To incorporate several magnons with the same momentum and spin, we can make a 
small modification to the integral equations. Multiple magnons give multiple copies of 
the densities p{x). Hence the integral equations now read 

(]r'n(r') 

Mj- / = -2 z M, ln(a; - A) + Cm, (3.26) 

J A X — X' 

where p{x) is the density for one magnon and 

E f \ 1 \ 
Cm = - + - 27m + 2Mz \n{A - A*). (3.27) 



2 \A-g A + g^ 

The solution is the same as in the single magnon case, provided we make the substitutions 

E J 6 , , 

M M ^ M ^ ^ 

The resulting shift of the dispersion relation is 

AS = -64sin^^cos(p-2^)^exp(-2^^/^f). (3.29) 
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This result is more transparent in the Q ^ hmit, where it simphfies to 



IQg sin^ f ^ n , ( JIM \ 

AS = ^ + 2.^) -p(-2i^)- 



A physical string would have p = 27r m/M, for some integer m. Hence we can remove 
the cos-factor by a suitable choce of n. This then reproduces AFZ's result for these 
special momenta, as well as the orbifold results of Astolfi et al. [23] . The physical strings 
with these M multiple magnons correspond to the limit of the GKP spinning string on 
5*5 [24] for M = 2, or their multicusp generalizations [25] when Q = 0. If Q 7^ 0, then 
the physical strings are the helical strings described in [34]. 

Note that the exponential suppression is not as large as for single magnons. Looking 



at (3.25), we see that the singularity as two magnons approach the same momentum is 



signifying this change in the suppression factor. 



4 Checking results — back to Heisenberg 

For Q ^ (? we expect the string result to approach the result from the one loop gauge 
theory calculation, i. e. from the Heisenberg model. For infinite spin J we have 

S-Q + —sm -. (4.1) 



From (3.21) the finite size correction in this limit is given by 



AS = -64 sin^ - cos(p - 20) ^ exp -2£^— (4.2 



^ cos(p - 20) ^ exp (^-2£ 

aP I nA\9^ f c,J + Q ■ 2P 



— 64sin - cos(p — 20)— exp I —2 — — — sin - ) , (4.3) 
2 Q \ Q 2 ^ 



where the phase is given by 



Q cot I _ ^ ^ ^ sinp 



E ^^-^-nn^iJ + Q) ^-Ttn. (4.4) 



It is instructive to reproduce (4.3) from the Heisenberg model. We start by quickly 
recalling the "Bethe string" solutions of the Bethe equations in the infinite J limit. 
Consider again the Bethe equation 

Uj-i/2) l}.Uj-Uk-i' 

Assume that Imuj > for some j. As J ^ 00 the left hand side will grow indefinitely. 
For (4.5) to be satisfied, there has to be a pole on the right hand side. We will have 
a pole if there is another root, Uj+i such that Uj — Wj+i —>■ i. If ImMj+i > then its 
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Bethe equation will also have a divergent left hand side. But its right hand side already 
has a zero since uj+i — uj —>■ —i. Hence to compensate, there must be a with 
Mj+i — i such that — > 0. If we continue finding roots in this way, 

eventually we will have a root Uk where \muk < 0. In this case, its Bethe equation will 
approach on the left hand side, so the right hand side will also have to be zero, hence 
we would need — Mfc+i — > i such that ^^-Uk-i+t _^ ^j^^ argument works in 

reverse from the Im Uj > case. 

Thus, as J — s> oo the Q Bethe roots will distribute themselves along a line parallel to 
the imaginary axis, each root separated by i. For the energy to be real, the configuration 
has to be symmetric under complex conjugation. The top root will then be Ui = Xq + 
i{Q — l)/2 and the bottom root uq = Xq — i{Q — l)/2. This distribution is shown in 
figure [6] 

To calculate the momentum p we use the equation 



e 



which we can solve for the parameter xq to get 



2 2 

The dispersion relation is given by 
Q 



Xo = — cot -. (4.7) 



S = Q+ig > — — — = Q+ig — — — - = 

Using the above solution for xq we get 



''\Uj + i/2 Uj-i/2) \ui+i/2 UQ-i/2) Q"^ + Ax, 

(4.8) 



£: = Q + |?sin^| (4.9) 

as expected. 

If we now consider J to be large but finite, the roots will not be exactly separated by 
i. To see how the distribution changes, suppose that ui is the top root in figure |6] and let 
Ui — uo = i + e. From the above argument we have that — as J — oo, hence we 

expect U2 — U3 = i + O(e^), — = i + C(e^), etc., until we reach a Uj that crosses the 
real axis. Once this happens the degree of separation reverses symmetrically. Hence the 
leading order finite size corrections will be governed by the small change in the positions 
of the top and bottom roots. This is similar to what we had for the finite gap equations, 
where the corrections corresponded to additional cuts at the ends of the condensates. 

Let us apply this to the case of a single set of bound magnons. Writing out the Bethe 
equation for ui = xq + i{Q — l)/2 we get 

x, + t{Q-l)) - e-l-2...(Q-2) " e ' ^ ' ^ 
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Figure 6: Bethe string 



where L = Q + J. For large L, the left hand side approaches exp (^ ^^^iQ ^ , and the 
parameter e can be calculated as 



Note that when comparing to (3.19), e ~ (5^6^*"^. Hence we expect nonvanishing finite 
size corrections already at leading order for the Heisenberg model. 

For finite L, xq will be shifted to x and the momentum will be shifted by 

^iAp ^ ^ + ^<3/2 + x-iQ/2 

x + iQ/2 x-tQ/2-ie*' ^ ' ' 

The shift in the dispersion relation can be calculated as 



AS = 2g^ — — + c.c. . (4.13) 

^ \x + iQ/2 + ie x + iQ/2 J ^ ' 

By adjusting x so that Ap = 0, we finally can express the correction to the dispersion 
relation as ^ 

A£ = -64sin^-cos(p-20)^exp ("-2^^^-^ sin^ ^"j . (4.14) 



This result from the Heisenberg model perfectly agrees with (4.3) in the Q ^ g limit of 
the finite gap solution. This procedure can be generalized to several sets of magnons, 
where one will continue to find agreement with the finite gap calculation. However in the 
case of condensates with density —in we need to be careful, since in the Bethe ansatz, 
it is not possible to have a Bethe string with density greater than —i [35] What we 
believe happens here is that there will be n nearby strings, with a small separation in 
the real direction between each string. 



^We thank K. Zarembo for comments on this point. 
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5 Results from TBA 



As a final check on our results, we consider the leading finite size corrections using the 
Thermodynamic Bethe ansatz . The TBA was first applied to = 4 by Ambjorn, Janik 
and Kristjansen [13], where they showed that it leads to corrections of order g^^ in the 
weak coupling expansion, precisely the order at which wrapping effects should appear 
[12]. More recently Janik and Lukowski applied the TBA to an explicit computation of 
the finite size corrections of a single fundamental magnon [36]. For another interesting 
application see [37]. 

Applying the TBA with more than one magnon would appear to be a hard problem. In 
order to find the complete leading order correction for M magnons, including prefactors, 
one would need to compute S-matrices for M+1 particles [36]. Although the final answer 



in (3.25) is quite simple, we have not yet attempted this computation. The S-matrices 
are fiavor dependent, but we expect this dependence to drop out when computing the 
leading order contribution. 

Instead we will focus on the more modest goal of reproducing our results for one 
dyonic giant magnon of charge Q and limit ourselves further by only finding the leading 
exponential term and not the prefactor |^ 

As was argued in [13,36], there can be two types of graphs that contribute to finite size 
effects, as was first shown for a relativistic field theory by Liischer [38]. These are called 
F- and /i-terms and are shown in figure [7j For the F-term graph the particle emits and 
absorbs an on-shell virtual particle that travels around the compact dimension. For the /x- 
term graph, the particle splits into two virtual particles that travel in opposite directions 
around the compact direction and recombine into the original particle. It turns out that 
the /x-term is responsible for the finite size corrections to the giant magnon energy [36] . 

The generalization of the /x-term graph in [36] is the creation/annihilation of two on- 
shell particles with charges n and Q', where n > 1 and we assume Q' 3> n. The charges 
are highest weights of SU{2\2) x SU{2\2), and so in order to conserve this nonabelian 
charge, they must satisfy Q — n<Q'<Q + n. Defining Q' = Q + n', we see that 
—n < n' < n. The exponentially suppressed correction then arises from the exp{ipcJ) 
factor that appears in the virtual particle propagator, where pc is the virtual particle 
momentum. The other virtual particle has momentum defined to be —6p, so on shell 
conservation of momentum gives Pc = P + Sp. Requiring energy conservation then leads 
to 



g2 + 16^2 sin2 E^hQ^ ^ ^Qg2 sin2 P^^J^2^ lQg2 sin2 ^-1 (5.1) 



To leading order (5.1) reduces to 



2 

which has the solution 



S^Sp"^ - An'Qg^ cot I5p + An^ = Q (5.2) 



O pnt - 9 I 

5p = 2n'^—^ + I ^Jn^Sl - n'2Q2 ^ot^ 2 (5.3) 



^As this paper was being prepared [28] appeared which does calculate the prefactor. 
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Figure 7: World-sheet diagrams corresponding to the F— and ^-terms. We have indicated the 
charge of the magnon as well as the virtual particle charges. 



We wish to minimize the imaginary part to lessen the exponential suppression. Clearly 
this occurs when n = 1 and n' = ±1 and so we find 



5p = ±2 



Qcotf 



£ 



(5.4) 



Both possibilities for n' will contribute to the /^-term graphs. If we now substitute (5.4) 



into exp{ipcJ), we obtain the exponential suppression in (3.21) as well as the phase in 
dslQl). 



6 Discussion 

In this article we have derived the leading finite size corrections for a general configuration 
of giant magnons on physical strings, that is when the total momentum is a multiple 
of 27C. We have found that the form of the corrections is particularly simple if all 
magnons have Q = 0. We have also shown that while the finite gap configuration for 
a general collection of giant magnons involves hyperelliptic functions, the leading order 
computation is found by reducing to a single cut and hence ordinary functions. 

One can go on and find the next order corrections as well. The method is straightfor- 
ward, albeit tedious. Again the computation can be reduced to single cut configurations 
and hence ordinary functions. One might also want to compare these results with other 
methods, for example by examining multi-soliton solutions in the sine-Gordon model. 

Another interesting question is to compute the quantization of the magnon momenta 
for large but finite J. Moreover, it would be interesting to see how a particular set of 
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magnons with quantized momenta connect onto string solutions with small values of J. 
Finally, it should be possible to compute the leading order corrections using the methods 
of Janik and Lukowski [36]. 
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